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Abstract 

It is argued that if cosmic rays penetrate into molecular clouds, the 
total energy they lose can exceed the energy from galactic supernovae 
shocks. It is shown that most likely galactic cosmic rays interacting 
with the surface layers of molecular clouds are efficiently reflected and 
do not penetrate into the cloud interior. Low-energy cosmic rays (E < 
1 GeV) that provide the primary ionization of the molecular cloud gas 
can be generated inside such clouds by multiple shocks arising due to 
supersonic turbulence. 


1 Introduction 

Gas phase chemical reactions in molecular clouds are predominantly ion- 
molecular reactions, and are catalyzed by H + ,H^, and ions, which 
are produced, in turn, by low-energy cosmic rays (E £ 1 GeV; see the re¬ 
view [1]). The primary ionization rate by cosmic rays required to sustain 
the fractional ionization and concentrations of these ions is estimated to 
be £ ~ ICG 1 ' s" 1 . A similar value is obtained for the upper limit to the 
primary ionization rate derived from the abundances of HD in molecular 
clouds, C A ICG 16 — ICG 17 s' 1 [2, 3]. It was first proposed in [4] that the 
fractional ionization in the cores of molecular clouds might exceed the value 
corresponding to ( ~ ICG 17 s _1 by a factor of five, which would require 
that ( be increased by a factor of 25. Based on the assumption that the 
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predominant magnetic field structure in molecular clouds is dipolar, it was 
shown in [5] that cosmic rays with energies E > 1 MeV/nucleon essentially 
freely penetrate into their interior, with ionization losses constituting about 
10% of their energy (see also [ 6 ]). For this reason, the above estimates of 
the primary ionization rate are considered to be valid for the interstellar 
medium as well. However, new data on the properties of molecular clouds 
and new theoretical analysis of their nature have recently cast doubts on the 
unconditionality of the conclusions of [5] and the general view that the pri¬ 
mary ionization rates by cosmic rays in molecular clouds and in the ISM are 
in agreement. We present here arguments supporting the idea that cosmic 
rays with energies E ^ 1 GeV are not able to penetrate into the interiors of 
molecular clouds, so that the presence of such particles there suggests that 
high-energy, nonrelativistic, nonthermal particles are generated within the 
molecular clouds themselves. This could occur, in particular, due to MHD 
shock waves arising in clouds during their formation from diffuse gas. 

Section 2 presents estimates of the total energy lost by galactic cosmic 
rays in molecular clouds, and shows that these energy losses exceed the total 
energy rate injected by supernovae into cosmic rays. Section 3 presents 
arguments based on the results of recent numerical modeling that suggest 
that the magnetic held structure in molecular clouds should differ from a 
simple dipolar held, and probably has a toroidal (helical) component in the 
plane perpendicular to the rotational axis of the cloud. In this case, most of 
the cosmic rays incident on the cloud should be reflected. Section 4 presents 
estimates of the efficiency of the acceleration of cosmic rays by MHD shock 
waves in molecular clouds. Finally, Section 5 summarizes our conclusions. 

2 Energy losses of cosmic rays in molecular clouds 

It can easily be shown that the dominant contribution to the primary ion¬ 
ization rate by galactic cosmic rays with an energy spectrum [7] 

3 x 10 4 

J(E) ~ --(cm 2 sr s MeV) -1 , E < 10 GeV (1) 

v ' (800 + E) 2 - 5 y w 

is provided by cosmic rays with energies E < 1 GeV (the energy E in (1) 
is in MeV); the spectrum has a break at energies E > 10 GeV, where its 
slope becomes steeper. The fraction of the total energy E = 7 m p c 2 lost by 
a cosmic ray particle (proton) as it passes through material with a hydrogen 
density along the path Nl is determined to order of magnitude as [ 6 ] 
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where A = ln[2?n e c 2 /3 2 //(l — /? 2 )] — (3 2 ~ 5 — 12 for energies from 1 MeV to 1 
GeV, A 722 = 10 _ 22 Afi, (3 = v/c, and c is the speed of light. As it penetrates 
into a cloud, a high-energy particle undergoes numerous scatterings, so that 
the total path length travelled by the particle in the cloud is L 4cR 2 /D, 
where R is the radius of the cloud and D is the diffusion coefficient for 
cosmic rays in the cloud. The total density of particles along the trajectory 
is Nl ~ 2NR/X N [ 6 ], where A is the characteristic correlation scale for 
turbulence in the magnetic field, A <C R, and N is the surface density of 
the molecular cloud. Observed molecular clouds display a roughly constant 
(independent of radius) surface density of IV = 1.5 x 10 22 cnr 2 (see the 
review [ 8 ]). Thus, 


5 e ~ (1.5 - 3) x l(r 4 ^^— , (3) 

Xfi = X/R; for = 100 MeV and Xr = 0.1, we find 5e ~ 0.01. 

The energy absorbed by all the molecular clouds in the Galaxy is 

OO 

EZc = S J SeEJ{E)cLE ~4x lO " 5 ^ 1 erg/s, (4) 

100 MeV 

where S is the total surface area of a molecular cloud. Adopting for the 
entire Galaxy the molecular cloud spectrum obtained in [9], 

— = AM' 1 ' 73 , (5) 

dM 1 1 1 

over the interval M = [M\ = 10 2 ,Af 2 = 10 6 ]Mq with the normalization 

factor A = 3.5 x 10 7 , we can find the total surface area of clouds 

m 2 

S = 4vr f R 2 ^-dM = 7 x 10 4f Vn“ 2/3 , (6) 

J dM ’ w 

Mi 

where n ~ 10 2 cm -3 is the mean gas density in the clouds, and a > 1 is a 
factor accounting the fact that the surface area of a cloud is greater than 
that of a sphere due to irregularities in the cloud boundary [10]. When 
calculating the normalization coefficient A, we assumed that the mass spec¬ 
trum of molecular clouds determined in [9] for the Perseus arm is valid for 
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the Galaxy as a whole, with the mass of molecular gas in the Galaxy being 
5 x 1O 9 M 0 . This yields 5~ 3 x 10 45 cr cm 2 , and E c £ IC — 2 x 10 42 <t/Aj? erg 
s^ 1 for the energy losses of cosmic rays in clouds. The production of energy 
in galactic cosmic rays in supernovae is E c g N ~ 10 i2 r/ui/ 30 erg s -1 , where 
r] < 1 is the fraction of the supernova energy that goes into cosmic rays (for 
which modern estimates yield i] ~ 0.1 [11]), and U\ / 3 o is the supernova rate 
in the Galaxy in units of explosions per 30 years. It can readily be verified 
that, for r) ~ 0.1 and \r ~ 0.1, the amount of energy in cosmic rays that 
is absorbed by molecular clouds is higher than the amount of cosmic-ray 
energy produced in supernovae. In fact, this conclusion can be stated more 
firmly: with the adopted assumptions, Eg r N <C E^ c , which implies that 
probably Ar <C 0.1. Indeed, if in reality the correlation length is deter¬ 
mined by viscosity, then A# {£/ -R) 3 / 4 , where £ is the mean free path of 
the particles (£ <C R ) [6]. Recent numerical modeling of MHD turbulence 
in molecular clouds [12, 13] suggests that the correlation length is 0.01-0.03 
of the scale on which energy is injected into the system, which in our case 
is the radius of the cloud. To resolve this inconsistency between the cosmic 
ray energy that could be absorbed in molecular clouds and the cosmic ray 
energy produced by supernovae, we must suppose that galactic cosmic rays 
interacting with molecular clouds are predominantly reflected, so that only 
a small fraction of the cosmic ray energy is lost. 

3 Reflection of cosmic rays from the surface layers 
of molecular clouds 

According to recent numerical modeling, molecular clouds represent tran¬ 
sient regions of enhanced density arising due to the action of converging 
flows of interstellar gas [14-17], whose lifetimes, tMC ~ 2 x 10' — 10 s yr, 
are probably determined by the activity of the stars born in them. In this 
picture, one of the main sources of turbulent energy in molecular clouds is 
the kinetic energy of the converging flows, although it is not ruled out that 
an appreciable amount of energy can be contributed by young stars born 
in the molecular cloud cores. During the formation of a molecular cloud - 
i.e., the transformation of atomic into molecular hydrogen in the process gas 
compression - some fraction of its rotational angular momentum is lost, but 
the remaining specific angular momentum of the cloud is still substantial: 
on average, the specific angular momentum in molecular clouds is a factor 
of four lower than the specific angular momentum of the diffuse interstellar 
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gas. This means that on the lifetime a cloud rotates about 30 times as it is 
compressed in the converging flows (we have adopted here a characteristic 
rotational velocity for the cloud 10 14 s 1 [18]). Thus, we expect the 
magnetic held in the outer parts of the molecular cloud to have a helical 
structure with an appreciable toroidal component (i.e., the component per¬ 
pendicular to the rotational axis of the cloud). This type of held structure is 
indeed observed in a number of cases (such as the cloud Lynds 1641) [19]. In 
this case, no matter what the angle at which a cosmic ray particle is incident 
on the molecular cloud surface, the pitch angle will reach its critical value 
as the particle penetrates into the cloud, causing the cosmic ray particle to 
experience mirroring. The characteristic time for the rehection of cosmic 
rays by this process is [5] 

'jm p fdlnBX ' 1 7 m p L B 

tm ~ 1 n ) ~ 

p \ OS J p 

where Lb is the length along the trajectory of the cosmic ray particle over 
which the regular component of the magnetic field varies, for which it is 
natural to adopt Lb ~ R. The characteristic diffusion time for the cosmic 
rays is 


R 2 

t-D ~ -jj- ( 8 ) 

If we suppose that this diffusion is determined by kinetic processes with 
cross section <r, ~ 3 x 10~ 18 cm 2 (the cross section for ionization losses), the 
diffusion time will be to = (IVcrj) 3 / 4 .R/c, so that 


to ~ {NaifR 


< 1. 


(9) 


If the diffusion is determined by resonance scattering on small scales of the 
order of the gyroradius of a cosmic ray proton, l ~ xr p [5, 6], the diffusion 
coefficient will be D ^ ci? 1 / 4 / 3 / 4 , and the ratio of the time scales becomes 


f / I \ 3 / 4 

r D ~ ( 5 ) « 1 (10) 

Thus, cosmic rays should be reflected by the enhanced magnetic field of 
a molecular cloud appreciably more rapidly than they can penetrate into 
the cloud interior. The fraction of cosmic rays that are able to penetrate 
into the cloud can be crudely estimated using the ratio t m /tD, which is to 
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order of magnitude 1CP 3 in both cases. This analysis forces us to conclude 
that the primary ionization rate of the molecular gas is substantially lower 
than the value that is required in order for the gas-chemical reactions in 
the cloud to be sustained: ( ~ 10 -1 ' s _1 . Indeed, galactic cosmic rays 
with the spectrum presented above provide to order of magnitude just this 
primary ionization rate. Even if about 10% of the high-energy particles can 
penetrate into the cloud, which could with some margin lead to agreement 
between E^c an d ^s"n- C would be at least an order of magnitude lower 
than the required value. This contradiction can be resolved assuming that 
turbulent motions in the molecular clouds are able to accelerate nonthermal 
particles to sufficiently high energies to sustain the required level of primary 
ionization. 

4 Acceleration of nonthermal particles in molecu¬ 
lar clouds 

Molecular clouds display well-developed turbulence, as a rule, supersonic, 
and apparently super-Alfvenic (see the discussions in [4, 12, 20]). It is, ac¬ 
cordingly, natural to expect the presence of numerous MHD shock waves in 
the cloud, as is observed in both numerical simulations (see the discussion 
in [21]) and in real molecular clouds [22, 23]. The characteristic velocities 
of MHD shocks in molecular clouds can reach ~ 10 km s -1 , in which case it 
is reasonable to expect that some fraction of their kinetic energy could be 
transformed into nonthermal particles accelerated in the shock fronts. The 
efficiency of particle acceleration under these conditions can be estimated 
using the theory developed in the studies [24-27], which are concerned with 
the acceleration of particles in the fronts of chaotic shocks in turbulent inter¬ 
planetary and interstellar fields. In this theory, the maximum energy that 
can be gained by particles is determined by the ratio of the time over which 
the acceleration region exists (the time for the shocks to travel through the 
distance separating them), t a ~ L s /u , and the time for the diffusion of the 
particles from the acceleration region, tp ~ L 2 S /D, where L s is the charac¬ 
teristic distance between the shock fronts and u is the shock velocity [26] 
(the accelerated particles leave the acceleration zone when t£>/t a < 1). The 
diffusion coefficient D = uA/3 is determined by the scattering of particles 
on small-scale inhomogeneities in the magnetic field with a characteristic 
transport free path A ~ C u L s {r p /L s ) 2 ~ u , where C u — 0.3, r p = cp/eB is the 
gyroradius of the particle, p is the momentum of the particle, and v is the 
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index of the spectrum of the small-scale fluctuations of the magnetic field 
(■ d{B 2 )/dk ~ k ", kL s 1 [26, 27]). We will estimate the maximum energy 
of the accelerated particles Em from the condition to/t a = 1. This yields 
for the spectral index v = 1.5 

Em ~ 5 x 10 3 GeV, (11) 

when Em < m p c 2 and 

Em ~ GeV, (12) 

when Em > m p c 2 ; here, = L s /R , i?i = R/(l pc), i?i = B/( 1 /jG). 
Thus, in clouds with characteristic parameters B ~ 10/uG and i? ~ 10 pc 
and with L s /i? ~ 0.01 — 0.03, we expect protons to be accelerated to ener¬ 
gies Em ~ 5 — 15 MeV. However, Em is very sensitive to the index of the 
magnetic Held fluctuation spectrum Em oc ( u/c j 1 /! 2- ") i n the relativistic 
limit and Em oc ( u/c) 2 E 2 ~ 1 ') j n the nonrelativistic case. Therefore, with 
a Kolmogorov turbulence spectrum, u ~ 1.7, and shock velocities u ~ 10 
km s _1 , the estimates of Em fall to values of the order of atomic ener¬ 
gies. From this point of view only detailed observations of the magnetic 
field fluctuation spectrum will enable confident conclusion of whether or not 
nonthermal particles can be accelerated to sufficiently high energies in molec¬ 
ular clouds. Recent numerical simulations of MHD turbulence in molecular 
clouds suggest that, in most cases, a Kolmogorov spectrum is established 
near wavenumbers kL/2ir 10, but the spectrum becomes appreciably non- 
Kolmogorov at longer wavelengths (exceeding 0.1 of the cloud radius), where 
it has an index near zero [12, 13]. Thus, the transition from large scales, 
corresponding to the energy injection scales, to the scales on which an in¬ 
ertial regime is established encompasses a broad wavenumber range. If we 
associate the possibility of generation of cosmic rays with long-wavelength 
components of the turbulent motions in the molecular cloud, the energy 
Em could reach values of 1 GeV or higher. However, both the steady state 
characteristic energy of the nonthermal particles E* and the spectrum of 
the generated cosmic rays are determined to a considerable extent by the 
injected particle distribution function - i.e., by those particles that, having 
entered in an active acceleration region, are able to undergo fairly frequent 
interactions with the shocks there - and by the fraction r/ of the injected 
particles in the flow as a whole that intersect the shock front. One source 
of particles that is usually considered is the particles that become extracted 
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from the thermal background during the scattering of fairly energetic pro¬ 
tons on forming collisionless shocks. This process is probably not efficient 
in molecular clouds, because the thermal energy of the particles here is too 
low. Winds from young stars and/or low-energy cosmic rays entrained in 
the converging flows forming the molecular cloud may serve as a source 
of such particles. In any case, the question remains open and requires a 
separate analysis; more detailed discussions in application to shocks from 
supernovae and to the diffuse interstellar and interplanetary media can be 
found in [24, 26, 28, 29]. The value of the fraction of interacting injected 
particles q is very uncertain. The injection rate at the bow shock in the mag¬ 
netosphere of the Earth is q ~ 10 -3 [30]. Estimates carried out in [31] for 
supernova remnants leave open a very wide range for q\ 10 -5 < i] < ICE 1 . 
If we adopt the lower limit for molecular clouds, the characteristic energy 
of the nonthermal particles proves to be too low. Based on the fact that 
the characteristic energy of the nonthermal particles should be such that 
their total energy W p ~ qrip 2 /2rn p does not exceed the kinetic energy of 
the generated turbulent motions Wt ~ pu 2 / 2, we arrive at the estimate 
E* 0.5 q~ l m p c 2 {u /c) 2 ~ 50 keV. The Larmor radius of protons with such 
energies is only r p ~ 3 x 10 9 cm in a magnetic field of B ~ WpG, so that 
the protons should primarily occupy a limited volume near the surface of 
the accelerating shock. 

To estimate the primary ionization rate provided throughout the molecu¬ 
lar cloud by the generated nonthermal particles, we suppose that only those 
particles that at times t ~ t a have diffusion lengths equal to half the distance 
between the shock fronts (i.e., ~ L S) r = 1/2) are able to occupy the cloud 
volume fairly uniformly and contribute to the ionization of the medium in 
the cloud. Adopting a power-law spectrum for the nonthermal particles, 
dNp/dE oc E~ q , and assuming that the total energy contained in these par¬ 
ticles is comparable to the energy of the turbulent motions, / EdN p ~ Wt , 
we obtain the primary ionization rate 

Em 

c= [ v{E)^a(E)dE ~ 3 x l(T 18 n s“\ (13) 

J dE 

Em /2 

for q = 2 and 


C ~ 3 x 10~ 19 n s” 1 , (14) 

for q = 2.5; here, n is the gas density in the molecular cloud. Thus, for 



a mean gas density n ~ 10 2 cm -3 [32], we find Q ~ 10 -1 ' — 10^ 16 s -1 , 
which is close to the value required to sustain the gas-phase chemistry in 
the molecular clouds. Note that larger values exceed the value expected 
from galactic cosmic rays. It is worth stressing in this connection that 
recent observations of some molecular clouds require primary ionization rates 
that cannot be provided by galactic cosmic rays. In particular, the value 
( = 5.6 x 10 -17 s _1 with an uncertainty factor of about three is preferred 
in [33]; i.e., the upper end of the allowed interval corresponds to £ ~ 10~ 16 
s _1 . However, even the most probable value exceeds the value expected 
from galactic cosmic rays by a factor of five. This can be considered as an 
additional argument for the need for additional sources of cosmic rays in 
molecular clouds and the possible generation of cosmic rays inside them. 

5 Conclusions 

We have shown the following: 

1. The energy lost by galactic cosmic rays in molecular clouds may 
exceed the total energy converted into cosmic rays by supernovae. 

2. An appreciable fraction of cosmic rays with energies E ~ 1 GeV can 
be reflected from the surface layers of molecular clouds through mirror effect. 
Thus, an additional source of nonthermal particles is required to sustain the 
gas-phase chemical processes occuring in these clouds. 

3. One such source is particles accelerated on the fronts of MHD shocks 

that arise during the formation of molecular clouds. The primary ionization 
rate that can be provided by such particles is £ 10“ 17 - 10 -16 s -1 . 
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